Quantum chaos in multicharged ions and statistical approach to 
the calculation of electron-ion resonant radiative recombination 



G. F. Gribakin, A. A. Gribakina, and V. V. Flambaum 

School of Physics, The University of New South Wales, Sydney 2052, Australia 

(February 2, 2008) 

Abstract 



We show that the spectrum and eigenstates of open-shell multicharged 
atomic ions near the ionization threshold are chaotic, as a result of extremely 
high level densities of multiply excited electron states (10 3 eV _1 in Au 24+ ) 
and strong configuration mixing. This complexity enables one to use sta- 
tistical methods to analyse the system. We examine the dependence of the 
orbital occupation numbers and single-particle energies on the excitation en- 
ergy of the system, and show that the occupation numbers are described by 
the Fermi-Dirac distribution, and temperature and chemical potential can be 
introduced. The Fermi-Dirac temperature is close to the temperature defined 
through the canonical distribution. Using a statistical approach we estimate 
the contribution of multielectron resonant states to the radiative capture of 
low-energy electrons by Au 25+ and demonstrate that this mechanism fully 
accounts for the 10 2 times enhancement of the recombination over the direct 
radiative recombination, in agreement with recent experimental observations. 
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I. INTRODUCTION 



In this paper we investigate the spectrum and eigenstates of a multicharged positive ion 
at energies close to its ionization threshold I. Using Au 24+ (I = 750 eV) as an example, 
we show that this spectrum is dominated by multiple electron excitations into a few low- 
lying unoccupied orbitals. As a result, it is extremely dense, with level spacings ~ 1 meV 
between the states of a given total angular momentum and parity J n . The electron Coulomb 
interaction induces strong mixing of the multiply-excited configurations, which leads to a 
statistical equilibrium in the system. The latter is similar to a thermal equilibrium, and 
variables such as temperature can be introduced to describe it. This enables one to use a 
statistical approach in the situation where a full dynamical quantum calculation is simply 
impossible because of the enormous size of the Hilbert space 10 5 for Au 24+ ). 

We apply this approach to the problem of radiative capture of low-energy electrons by 
multicharged positive ions, and show that in these systems the contribution of resonant mul- 
tielectronic recombination that proceeds via electron capture into the multiply- excited com- 
pound states, is responsible for high recombination rates, much greater than those expected 
from the simple direct radiative recombination. Our calculation resolves quantitatively the 
long-standing puzzle of huge enhancements of the electron-ion recombination rates, and es- 
sentially removes the "enormous discrepancies between theoretical and experimental rate 
coefficients" (Hoffknecht et al. 1998). The situation here turns out to be similar to the 
radiative neutron capture by complex nuclei [(n, 7) reaction] where the resonance mecha- 
nism involving the compound nucleus states is also much stronger than the direct capture 
(Flambaum and Sushkov 1984, 1985). 

So far the enhancement of the recombination rates at low electron energies ^ 1 eV has 
been observed for a number of ionsQ. Its magnitude ranges from a factor of about ten for 
Ar 13+ (Gao et al. 1995), Au 50+ and Pb 53+ (Uwira et al. 1997a), and U 28+ (Miiller and Wolf 
1997), to over a hundred for Au 25+ (Hoffknecht et al. 1998). This enhancement is sensitive 
to the electronic structure of the target, e.g., the recombination rates of Au 49+ and Au 51+ are 
much smaller than that of Au 50+ (Uwira et al. 1997a). For few-electron ions, e.g., C 4+ , Ne 7+ 
and Ar 15+ the observed rates are described well by the sum of the direct and dielectronic 
recombination rates (Schennach et al. 1994, Zong et al. 1997, Schuch et al. 1997). In more 
complicated cases, like U 28+ or Au 25+ , the questions of what are the particular resonances 
just above the threshold and how they contribute to the recombination "remain a mystery" 
(Mitnik et al. 1998). 



Apart from the enhancement at eV energies due to many-electron processes, which is the subject 
of our work, there is another specific enhancement at electron energies below 1 meV. This enhance- 
ment increases with the charge of the ion, and is observed for all ions including fully stripped ones, 
see Gao et al. 1997 and Uwira et al. 19976, and we do not consider it here. 
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II. SPECTRUM AND EIGENSTATES OF AU 24+ 



Let us consider the problem of electron recombination on Au 25+ . Due to electron cor- 
relations the low-energy electron can be captured into an excited state of the compound 
Au 24+ ion. This system is the main object of our analysis. Au 24+ has 55 electrons. Its 
ground state belongs to the Is 2 . . . 4/ 9 configuration. Figure [l] shows the energies of its 
relativistic orbitals nlj obtained in the relativistic Hartree-Fock calculation. All orbitals 
below the Fermi level, Is to 4/, were obtained in the self-consistent calculation of the Au 24+ 
ground state. Each of the excited orbitals above the Fermi level - 5s, bp, etc., was calcu- 
lated by placing one electron into it, in the field of the frozen Is 2 . . . 4/ 8 core. The energy 
of the highest orbital occupied in the ground state is £4/ 7/2 = —27.9 a.u. This value gives 
an estimate of the ionization potential of Au 24+ : I ~ | £ 4/ 7 / 2 l- Our relativistic configuration- 
interaction (CI) calculation of the ground states of Au 24+ 4/ 9 and Au 24+ 4/ 8 shows that they 
are characterized by J = y and 6, and their total energies are —18792.36 and —18764.80 
a.u., respectively. Thus, the ionization threshold of Au 24+ is / = 27.56 a.u.= 750 eV, in 
agreement with Hoffknecht et al. 1998. 

The excited states of the ion are generated by transferring one, two, three, etc. electrons 
from the ground state into empty orbitals above the Fermi level (Fig. [1]), or into the 
partially occupied 4/ orbitals. We are interested in the excitation spectrum of Au 24+ near 
its ionization threshold. This energy (27.5 a.u.) is sufficient to push up a few of the nine 4/ 
electrons, and even excite one or two electrons from the Ad orbital. However, the preceding 
4p orbital is already deep enough to be considered inactive. Thus, we treat Au 24+ as a 
system of n = 19 electrons above the frozen Kr-like Is 2 . . . 4p 6 core. Note also that infinite 
Rydberg series corresponding to the excitations of one electron in the field of Au 25+ belong 
to the single-particle aspect of the Au 25+ + e~ problem, and we do not consider them here. 

The number of multielectron states obtained by distributing 19 electrons over 31 rela- 
tivistic orbitals, Ad^/2 through to 7g$/2, is enormous, even if we are only interested in the 
excitation energies below 27.5 a.u. It is impossible to perform any CI calculation for them. 
However, there is another simpler way to analyse the spectrum. The scale of the configu- 
ration interaction strength is determined by the two-body Coulomb matrix elements which 
transfer electrons between different configurations. Their typical size in neutral atoms is 
~ 1 eV, and in Au 24+ it is about 1 a.u., which is roughly 25 times greater, due to the 
smaller radius of the ion. This scale is much smaller than I. Configuration mixing aside, 
the CI does not shift the mean energies of the configurations. Therefore, we can construct 
the excitation spectrum of Au 24+ by calculating the mean energies Ei of the configurations, 
and the numbers of many-electron states Nj. within each of them: 

E t = E com + £ e a n a + £ Wfl( "» ~ , (1) 

a<b 1 + d <* 

Ni = U u 9a - (2) 
a nj{g a -n a )\ 

where n a are the integer orbital occupation numbers of the relativistic orbitals in a given 
configuration (J2a n a — n), e a = (a\H corc \a) is the single-particle energy of the orbital a in 
the field of the core, g a = 2j a + 1, and U a b are the average Coulomb matrix elements for the 
electrons in orbitals a and b (direct minus exchange): 
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Here R abba is the two-body radial Coulomb integral of A multipole, and 6 P — 1 when Z a + Z& + A 



is even, and otherwise. The mean energy of the lowest configuration 4d 10 4ff, 2 4fj , 2 obtained 



from Eq. ([]]) is just 0.28 a.u. above the CI ground state. 

Using Eqs. dH)-(0) we find that there are 9000 configurations within 35 a.u. of the Au 24+ 
ground state. They comprise a total of 2.1 x 10 8 many-electron states. If we allow for about 
10 different values of J, 2 J + 1 values of J z and the two parities, there would still be about 
5 x 10 5 states in each J n manifold. In Fig. |^ we show the total density of states for Au 24+ 
as a function of y/~E, where E is the excitation energy of the system above the ground state. 
It is obtained by smoothing out the small-scale fluctuations of the level density 

p(B)= y £N i 5(E-E i ) (4) 

by folding it with a Gaussian with 1 a.u. variance. In reality this averaging is done by 
the interaction and mixing of the configurations (Flambaum et al. 1994), but the result is 
expected to be the same. The inset on Fig. ^] presents a break-up of the total density near 
the ionization threshold into the densities of states with given J: p(E) = ^2j(2J + l)pj(E). 
The most abundant values are J = | to ^. For a given parity the density of such states at 
E « I is pjn pa 3.5 x 10 4 a.u., which corresponds to the mean level spacing D = l/pj* ~ 1 
meV. Figure |2] demonstrates the characteristic p oc exp(av^E) behaviour of the level density 
predicted by the Fermi-gas model (Bohr and Mottelson 1969), where a is related to the single- 
particle level density at the Fermi level g(sp) as a = [2ir 2 g(eF )/3] 1/ ' 2 . g{eF) = 3a 2 /2n 2 . We 
obtain an accurate fit of the level density at E > 1 a.u. by using a Fermi-gas model anzats 

p(E) = AE~ u exp(aVE) , (5) 

with A = 31.6, v = 1.56, and a = 3.35. The corresponding value of g{ep) = 3a 2 /27r 2 = 1.7 
a.u. is close to what one obtains from the Hartree-Fock orbital spectrum in Fig. p]. The other 
two parameters are different from the non-interacting Fermi-gas model values A = 1/ -\/48 
and v — 1. The latter values in fact lead to strong underestimation of the level density. For 
most abundant J n states the density pjn(E) is given by Eq. (|5[) with Aj* ph 0.15. 

At first sight the huge level density makes the spectrum of Au 24+ enormously compli- 
cated. On the other hand, this complexity enables one to analyse the system using statistical 
methods. The interaction between multiply-excited configuration states mixes them com- 
pletely, and they loose their individual features. In this regime the spectral statistics become 
close to those of a random matrix ensemble, the eigenstates cannot be characterized by any 
quantum numbers except the exact ones (energy and J 7 "), and the orbital occupation num- 
bers deviate prominently from integers. This regime can be described as many-body quantum 
chaos. We have extensively studied it in direct numerical calculations for the rare-earth atom 
of Ce - a system with four valence electrons (Flambaum et al. 1994, 1996, 1998a, 19986, 
Gribakina et al. 1995). 

The strength of the configuration mixing is characterized by the spreading width r spr . 
For a configuration basis state $fc with energy E^ it defines the energy range [E — E^l ^ r spr 
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of eigenstates in which this basis state noticeably participates. By the same token it shows 
that a particular eigenstate \I/ = J2kC k &k contains a large number N ~ T spr /D of principal 
components - basis states characterized by C k ~ 1/ VN. Outside the spreading width C k 
decrease. This effect is usually referred to as localization. Apart from this, C k behave 
closely to Gaussian random variables (Flambaum et al. 1994). The effect of spreading is 
approximated well by the Breit-Wigner shape of the mean-squared components (Bohr and 
Mottelson 1969) 

i r 2 /4 

The normalization Ylk C\ = 1 yields iV = 7iT /2D. In systems with small level spacings D the 
number of principal components N can be very large. It reaches several hundreds in Ce, and 
can be as large as 10 6 in complex nuclei. At \E k — E\ > T Eq. (|6|) gives Cj;(E) oc l/(E k — E) 2 , 
which corresponds to the simple first-order perturbation theory dependence with constant 
mean-squared mixing matrix elements. In real systems the mixing between distant (in the 
sense of their unperturbed energies) basis states is usually suppressed. Accordingly, the 
Hamiltonian matrix is characterized by certain handedness, i.e., the off-diagonal matrix 
elements iiy decrease as one moves away from the main diagonal i = j (Gribakina et al. 
1995). This causes a faster, close to exponential, decrease of the mean-squared components 
at large \E k — E\ (Flambaum et al. 1994). 

In Fig. |H we illustrate the behaviour of the eigenstate components by the results of a 
CI calculation which includes just two odd configurations of Au 24+ with energies close to 
the ionization threshold: 4/|/ 2 4/| / 2 5pi/2 5^3/2 5/V /2 and 4f^ 2 ^f7/2^Pi/^^3/^97/2- These two 
configurations produce a total of 143360 many-electron states with J ranging from | to 
17.5. As an example we present the results obtained by diagonalization of the Hamiltonian 
matrix for J n = y . This total angular momentum value is among the most abundant 

in the spectrum, as there are 1254 J n = y states. The mixing of the two configurations 
included is practically complete, since the weight of each configuration in every eigenstate 
is close to 50%, Fig. |1 Shown in the upper part of Fig. ^ are the components of the 590th 
eigenstate from the middle of the two-configuration spectrum. Both the fluctuations of C k 
as function of basis state k, and the localization of the eigenstate components in the vicinity 
of the corresponding eigenvalue (E = 27.51 a.u. above the Au 24+ ground state) are evident. 

A Breit-Wigner fit of the mean-squared components yields N = 975 and r spr = 0.50 
a.u., see lower part of Fig. |3|. When the calculations are performed for one of the above 
configurations, N is about two times smaller, but r spr is practically the same. The spreading 
width is related to the mean-squared off-diagonal Hamiltonian matrix element and the mean 
level spacing as r spr ~ (Bohr and Mottelson 1969). It is known to be a robust 

characteristic of the system. When more configurations are included, both D and Hfj 
decrease, whereas r spr does not change much. If one could do a full-scale CI calculation near 
the ionization threshold of Au 24+ one would obtain eigenstates with iV = (7r/2)r spr pj7r ~ 
3 x 10 4 principal components. 
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III. STATISTICAL APPROACH 

The spreading of the basis states due to configuration interaction introduces natural 
statistical averaging in the system. Based on this averaging, a statistical theory of finite 
Fermi systems of interacting particles can be developed (Flambaum and Izrailev 1997a, 
19976). It enables one to calculate various properties of the system as sums over the basis 
states, without actually diagonalizing the Hamiltonian matrix. For example, the mean 
orbital occupations numbers can be obtained as 

n a (E)=Y,Cl(EH k) (7) 

k 

where n^> is the occupation number of the orbital a in the basis state k. To demonstrate 
how it works we have used a simple Gaussian model spreading 

Cl(E) oc exp 

and calculated the mean orbital occupation numbers as functions of the excitation energy 
E using = 1 a.u., Fig. j|. Of course, in our calculation we sum over the configurations, 
rather than the actual many-electron basis states, and use their mean energies and weights 
given by Eqs. (|) and (|), cf. Eq. (g). 

The oscillatory dependence with the period of about 3-4 a.u. is due to the shell structure 
of the Au 24+ ion, Fig. [I]. As the excitation energy increases the oscillations die out. Apart 
from this the occupation numbers of the orbitals below the Fermi level (4d and 4/) decrease, 
and those above it (5s, 5p, etc.) increase, as one would expect in a Fermi system. It seems 
very natural to try to describe this behaviour in the spirit of statistical mechanics, by intro- 
ducing temperature and applying the standard Fermi-Dirac (FD) distribution (Flambaum 
et al. 19986). Temperature has long been used to describe highly excited nuclei, and the 
question of thermalization was investigated recently in numerical calculations for the s — d 
shell nuclear model (Horoi et al. 1995, Zelevinsky et al. 1996). Of course, temperature can 
always be used to describe the equilibrium of a macroscopic system that contains a large 
number of particles, or to describe a small system interacting with a heat bath. In what 
follows we are going to see if the notion of temperature can be applied to our isolated system 
with a small number of active particles. The total number of electrons in Au 24+ is quite 
large, however most of them are inactive at the excitation energies at or below the ionization 
threshold. 

The formula for the single-particle occupation numbers v a = n a /g a (0 < v a < 1) 

1 

1 + exp[(e - /i)/T] 

at a given temperature T and chemical potential /i depends on the single-particle orbital 
energies e a . These energies are well defined for non-interacting particles in a given potential. 
For interacting particles (electrons in an atom or ion) one can introduce single-particle 
orbitals and energies using a mean field approximation, e.g. the Hartree-Fock method. 
From this points of view we could use the orbital energies e^ F found in the mean field of 
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the Au 24+ ground state Is 2 . . . 4<i 10 4/ 9 , Fig. |l|. However, they may only be suitable at low 
excitation energies, when the mean field is close to that of the ground-state Au 24+ . 

As the excitation energy increases the orbital occupation numbers change noticeably, as 
shown by Fig. [|. This gives rise to a change of the mean field, and as a result, the orbital 
energies are shifted by 

5e a (E)=Y,U ab 5n b (E) , (10) 

b 

where 5n b = nb(E)—nb(0) is the difference between the occupation numbers at energy E and 
in the ground state at E = 0. Using our numerical energy- dependent occupation numbers 
we find the energy dependence of the orbital energies, shown for a few low-lying orbitals in 
Fig. g. With the increase of the excitation energy the electrons are transferred into higher 
orbitals which have larger radii. Accordingly, the electron cloud becomes more diffuse, the 
screening of the nuclear potential is reduced, and the orbital energies go down. This effect 
is especially strong for the inner Ad and Af orbitals. As we will see below the shift of the 
lower orbital energies is comparable with the temperature of the system, and it has to be 
taken into account when applying the FD formula. 

In Fig. [7] we present the single-particle occupation numbers at five different excitation 
energies, as functions of the shifted orbital energies 

e a {E) = e a (0) + 6e a (E) , (11) 

where we take e a (0) = e^ F . The numerical values agree well with the FD distribution Eq. 
([]), obtained by the least-square fits of the temperature T and chemical potential p (solid 
lines). Figures |8| and |^ present the dependence of p and the "Fermi-Dirac temperature" T 
on the energy of the system (solid circles). 

Occupation numbers aside, the relation between the temperature and energy can be 
defined by the level density p(E), Eq. (|), through the canonical average 

Fm _ Je- E / T Ep(E)dE _ Y.jEjN^rr 
1 } Je- E / T p{E)dE EiNie- E >/ T ' 1 ] 

or from the statistical physics formula 

1 d\n\p(E)} . . 

T =^P< < 13) 

using the smooth fit (||). The latter yields T ~ 2y/~E/a, or E oc T 2 , characteristic of the 
Fermi systems. Figure [| shows that for the energies above 3 a.u. all three definitions of 
temperature give close values. As is known, the expansion of the chemical potential in 
a Fermi system at small temperatures starts with a T 2 term (Landau and Lifshitz 1969). 
Accordingly, its shift from the ground-state value is proportional to the energy. Indeed, a 
simple linear fit p = —27.6 — 0.094.E closely follows the numerical values in Fig. [8]. 

If we use T(E) given by the canonical definition and the linear fit of p, together with 
the orbital energies [Eq. ([TTj)] , the FD formula gives smooth energy dependencies of the 
occupation numbers, see Fig. [5[ They reproduce the behaviour of the numerical occupation 
numbers averaged over the shell-structure fluctuations. 



7 



IV. DIRECT AND RESONANT RECOMBINATION 



Let us now estimate the direct and resonant contributions to the recombination rate 
of Au 25+ . The direct radiative recombination cross section is estimated by introducing an 
effective ionic charge into the Kramers formula, which describes radiative electron capture 
in the Coulomb potential, see e.g. Sobelman 1992, 

ff(d ) = (^ 2R y d ) 2 fl4 ) 

n 3^3c 3 ne(Z?Ryd + n*e) ' 



where e is the initial electron energy, n is the principal quantum number of the final elec- 
tron state, and atomic units are used (Ryd = | a.u.). If we are interested in the total 
recombination cross section the sum over n must be calculated, 

^ (d) = E-i d) - (15) 

n 

Due to the n~ l factor in Eq. (|14]) this sum diverges logarithmically, until values of n ~ ra m ax 
are reached, where n^ ax e = Z 2 Ryd, after which it converges rapidly. With the logarithmic 
accuracy the result is given by 

oto^J^Wp jKf^ f (16) 
3^ e 1 \n \ e j ' K ' 

where no is the principal quantum number of the lowest unoccupied ionic orbital, which 
determines the lower limit in the summation over n. Using Zi = 25, uq = 5 for electron 
recombination with Au 25+ , and choosing a small electron energy of e = 0.1 eV we obtain 
0-( d ) k 7 x 10~ 17 cm 2 . This corresponds to the recombination rate of A = av = 1.3 x 10~ 9 
cm 3 s _1 , which is two orders of magnitude smaller than the experimental A = 1.8 x 10~ 7 
cm 3 s _1 at this energy (Hoffknecht et al. 1998). 

The electron energy of 0.1 eV is equal to the transversal temperature of the electron 
beam in the experiment, whereas the longitudinal temperature is much smaller, 1 meV. 
Therefore, to make estimates of the recombination rates at this and higher energies one can 
use the the cross sections without averaging over the Maxwellian velocity distribution. It 
is also important that the energy dependence of the experimental recombination rate is in 
agreement with that of the direct radiative capture for electron energies 1 meV< e ^ 1 



eV. The latter is basically given by the 1/e factor in Eq. (16). The experimental data of 
Hoffknecht et al. (1998) is reproduced well by the direct radiative rate multiplied by a factor 
of 150. 

The cross section of resonant radiative capture averaged over the resonances is (Landau 
and Lifshitz 1977) 

where T 7 and r e are the mean radiative and autoionization (or elastic) widths of the reso- 
nances, D is the mean resonance spacing, and we drop the statistical weights of the initial 
and intermediate ionic states. The relation T e ^> T 7 is usually valid for a few lower partial 
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waves, where the electron interaction is stronger than the electromagnetic one. Equation 
(|17|) is written for the electron s-wave, and the contributions of higher electron partial wave 
contain an extra factor (21 + 1). 

The radiative width of the resonant state at energy E »s I is found by summing the 
partial widths for all lower- lying states E' = E — u, 

r 3 /-'4cu 3 K| 2 

r 7 « L oLa - w)dw , (18) 



2 J + 1 7o 3c 3 

where the factor 3 accounts for J' = J, J ± 1, and d w is the reduced dipole matrix element 
between the many-electron states. Because of the chaotic structure of these states d u is sup- 
pressed compared to the typical single-particle matrix element d : d u ~ d /\/N (Flambaum 
and Sushkov 1984a, Flambaum et al. 1994, 1996). This estimate for systems with dense 
chaotic spectra in fact follows from the dipole sum rule: the number of lines in the spectrum 
is large, oc D^ 1 oc N, consequently, the line strengths are small, |eL| 2 ~ \do\ 2 N _1 . 

The integrand in Eq. (|i~8|) peaks strongly because of the competition between the u 3 
factor, and the level density p (I — to) that drops quickly as we go down from the threshold, 
see Eq. (|5|). As a result, the integral can be evaluated by the saddle-point method. Using 
the statistical estimate of d u we obtain 



(2J + i) c %r 



spr 



where ujq = G\/l/a corresponds to the maximum of the decay photon spectrum in Eq. 
This cross section has the same energy dependence as a^ d \ Hence, it is also in agreement 
with the energy dependence observed in the experiment, and we can estimate its magnitude 
at one particular electron energy, e.g., 0.1 eV. To do this we use a simple estimate of the 
single-particle dipole matrix elements in the ion with charge Zf. d ~ Zf , together with 
2 J + 1 w 10, and substitute r spr = 0.5, u = 9.4, and pj^[l — u ) = 2.5 x 10 3 a.u. into 
Eq. (P|). At e = 0.1 eV this gives a (r) = 7 x 10" 16 cm 2 , therefore, a® / a® = 10, and we 



obtain a factor of ten enhancement over the direct recombination due to radiative capture 
into multiply excited resonant states (the corresponding radiative width is T 7 = 2 x 10~ 7 
a.u.). It comes from the large effective number of final states in the radiative width in Eq. 
(|i"8|) (numerically T 7 ~ 2 x 10~ 7 a.u.). This enhancement has been obtained for the electron 
s-wave. The contributions of higher electron partial waves are similar to Eq. (|lj]) times 
(21 + 1). Therefore, a few lower partial waves (s, p, d) produce resonant cross section values 
10 2 times greater than a^ d \ which matches the experimentally observed values. With the 
increase of the orbital angular momentum I of the electron the capture width T e becomes 
smaller than the radiative width, and the contribution of the higher partial waves to the 
resonant cross section is suppressed. 



V. CONCLUSIONS 

In summary, the resonant radiative capture mechanism fully explains the strongly en- 
hanced recombination rates observed for eV electrons on multicharged ions. Its origin is in 
the high level densities of chaotic multiply-excited electron states in multicharged ions. The 
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size of the enhancement is sensitive to the electron structure of the ion, which determines 
the level density. We have shown that a statistical approach can be applied to the analysis 
of this complex system. One can also use a statistical theory to calculate mean-squared ma- 
trix elements between multiply excited chaotic states in terms of single-particle amplitudes, 
occupation numbers, r spr and D (Flambaum and Vorov 1993, Flambaum et al. 1994, 1996), 
and obtain accurate quantitative information about the processes involving chaotic states 
and resonances. At higher electron energies the resonant capture proceeds via so-called 
doorway states (Bohr and Mottelson 1969) - simple dielectronic autoionizing states, which 
are then "fragmented" into the dense spectrum of multiply-excited resonances (see Mitnik 
et al. 1998, Flambaum et al. 1996 and Refs. therein). 
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FIG. 1. Electron orbitals of Au + from the relativistic Hartree-Fock calculation. 
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FIG. 2. Level density in Au . Full circles connected by dotted line to guide the eye is the 
result of our numerical calculation. Solid line is the analytical fit, Eq. (|5|). The inset shows the 
densities of states with different J near the ionization threshold E = I. 
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FIG. 4. Weights of the ^fh 2 ^fy2^Pi/2^>Pz/2^f7/2 configuration in the J n = ^ eigenstates 
obtained in the two-configuration calculation. Note that for a few lower eigenstates the weights of 
this configuration are small, because its mean energy is about 0.03 a.u. higher than that of the 
other configuration. 
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FIG. 5. Energy dependence of the single-particle occupation numbers n a /g a calculated numer- 
ically from Eq. ([?]) for Ad, 4/, 5s and bp orbitals: solid and open circles correspond to j = I ± ^ 
subshells, respectively. Solid and dashed lines (j = I ± ^> respectively) show the results obtained 
from the FD formula using the energy-dependent orbital energies and chemical potential, and the 
canonical relation between the temperature and the excitation energy (solid line in Fig. |9|). 
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FIG. 6. Shifts of the single-particle orbital energies, Eq. ([[(]), as functions of the excitation 
energy for the 44/2,5/2 (long dash), 4f 5/2 j/ 2 (solid), 5s (dotted), 5p 1/2j3 /2 (chain), and 557/2,9/2 
(short dash) orbitals. 
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FIG. 7. Orbital occupation numbers in Au 24+ calculated numerically from Eq. (0) at excitation 
energies E = 1, 4.5, 9.5, 17 and 27.5 a.u. (solid circles), and the Fermi-Dirac distributions (solid 
line) with temperature T and chemical potential /i chosen to give best fits of the numerical data. 
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FIG. 8. Chemical potential obtained from the FD distribution fits of the occupation numbers, 
Fig. [7|, as a function of the excitation energy of Au 24+ (solid circles). Solid line is a simple linear 
fit /J, = -27.6 - 0.094E a.u. 
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FIG. 9. Temperature vs energy for Au 24+ . Solid line - canonical definition, Eq. (12); dashed line 
- statistical physics definition, Eq. (|l3|), which uses the density fit (||); solid circles - Fermi-Dirac 
fits of the occupation numbers. 
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